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1. Introduction 


The purpose of this paper is to study the distribution in space and energy 
of neutrons, slowed down by elastic collisions with the atomic nuclei of a 
moderator assumed to be homogeneous and infinitely large. The classical method 
of tackling this problem is to apply to the transport equation a Laplace trans- 
formation with respect to the logarithmic energy variable and a Fourier trans- 
formation with respect to the spatial coordinates. In the stationary and spher- 
ically symmetrical case, when in particular the source function is supposed to 
be a function only of the distance R from the centre of symmetry and of the 
energy HH, one obtains for the number of collisions per unit volume, per unit 
energy range Fy (R, EF) 

1 ny t+100 +00 
(1) F(R, B)-— eq | due | dowe®y(o, 1). 


m— toe Fa 


Here the logarithmic energy variable x is defined as x= log rae where Ey is 


a constant. The path of integration in the 7-plane is a straight line to the 
right of the singularities of the integrand. 

Equation (1) is the starting point of our investigation. The general theory 
leading to this equation was first published by WaLuER [8], [9]. Similar methods 
were also worked out in U.S.A. during the war, and appear in the com- 
prehensive review by Marsuaxk[6]. For details of the general theory we refer 
to these papers as well as a paper by VerpE and WricKx[10]. We only point 
out here the underlying assumptions made in this theory which thus must be 
fulfilled in order that our basic equation (1) be valid. Inelastic scattering is 
neglected; the elastic scattering is assumed to be isotropic in the centre of 
gravity system of the neutron and the atom; the atoms of the medium are 
supposed to be at rest before a collision with a neutron occurs and their binding 
to other atoms are neglected. Besides these assumptions and those mentioned 
before, we also assume for simplicity that the neutrons are emitted from a point 
source, with an initial energy equal to Hy. Further we neglect absorption, though 
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the way of treating a more general case without this restriction will also be 
shortly discussed. 

In general, the computation of F(R, Z) consists of two problems; firstly the 
calculation of the function fy (@, 7) in (1) and secondly the evaluation of the 
integral. However, in the general case when the mean free path A of the | 
neutrons is supposed to be a function of the energy # it is difficult to calculate 
fo(w, y) for all values of the variables necessary for the evaluation of (1) for 
any value of, for example, R. For small values of R (and of £) Fo (R, £) 
can be easily computed even in the general case; this formula of Fy (R, #) con- 
stitutes the well-known age approximation.! However, we intend to calculate 
the neutron density for any value of R, with the only restriction on the values 
of the parameters being imposed on £ (or x), and, for simplicity, we shall there- 
fore in this paper consider only the case when A is constant. Then, for every 
moderator R is measured with A as unity of length. The assumption of a con- 
stant A within the slowing-down energy interval is fulfilled with good approxima- 
tion only for sufficiently low values of Ey, roughly of the order 0.1 Mey. For 
higher Ey we may try to take some mean value for 4. If this mean value is 
suitably chosen, the error introduced by this approximation is considerably 
reduced and we may expect a qualitative agreement, at least. 

It is then easy to get an expression for fy(, 7), and we shall use that 
derived by WaLuER [8], [9], having the form of a continued fraction 


(2) fol, ") = 5 Gt atTB : 
hy (n) + 4@?/15 
ha (q)+9@?/35 
his (iq) 
with 
(9) ae on 


where @ is the number of neutrons emitted per unit time [9]. 
Further 


(4) hn (n) =1—en (n) (n=0, il. 2, aie “Je 


_ The expression of ¢n(y) has also been derived by Wa.rer [8], [9]. Making a 
simple transformation in his formula we get? 


M+1) £ G 2 : Y jee. 
(5) én (n) seul): | Px (11) il (4) ] ae 5=) 0; for M L 
| —lsdor Wie I 
with 


* See for instance MarsHak [6]. 
* Cf also Wicx [11]. 
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~, wt (uw? + M?—1)3 

(6) G (u) Mal 

Here MW is the mass of the scattering atomic nuclei measured in units of the 
neutron mass, P» (i) the n:th Legendre polynomial. Thus we have, for simplicity, 
assumed the moderator to contain one single slowing down element, the gene- 
ralization to the case of a mixture is, however, very easily done [9]. Finally 
from the value of the collision number Fy)(R, £) we get the simple expression 
of the neutron density No(R, H)=v~! F(R, EZ) (v= neutron velocity). 

Before we continue with this investigation we shall briefly sum up some 
earlier calculations of the neutron density beyond the range of validity of the age 
approximation. WALLER[9]| was the first to make such a calculation based on the 
general theory referred to above. He introduced the saddle point method im 
this problem. MarsHaxk[6] proposed a method of successive approximations, 
giving a series development for the neutron density. This series, however, 
diverges at distances essentially beyond the range of validity of the age theory. 
PuaczeEK! has made a partial summation of this series in order to use it at 
larger distances. However, the term he obtained is contained in WALLER’s 
formula. In a paper by the present author [4] WaLiEeR’s method was used 
and extended to higher approximations. Important progress has been made by 
Wick [11]. He has used the saddle-point method and derived the asymptotic 
formula of the neutron density for large distances. He has also described a 
numerical method of calculating the neutron density at any distance. In this 
paper we shall develop a method making it possible, by means of analytic con- 
tinuation, to obtain analytic formulas of the neutron density at all distances. 


2. Performance of one integration in (1) 


We first summarize some facts, well-known from the papers cited above, 
concerning the function fy(@, 7). In an w-plane cut along the imaginary axis 
from —ico to —? and from +7 to +700 fy is a one-valued analytic function 
of w with no other singularities than isolated poles provided 7 is chosen so 
that the relation 


(7) lim ¢p (4) =0, i.e. lim hy (y)=1, 


N+? n+ on 


holds. This is always valid in the right halfplane of 7. In the -plane fp has 
poles for certain values of w, occurring, of course, in pairs with opposite signs 
and being in general functions of 7. These poles are intimately connected with 
the zeroes of the functions h»(y) so that small values of the poles w always 
correspond to 7-values near to such a zero. One of the two poles corresponding 
to the zero with the largest real part, the root 7=1, of the equation ho (7) =0, 
gives the largest contribution to the integral (1) and is called the “main” 
pole ap. 

Now, for R{n}—co? which implies hy(n)>1, fo has a very simple limiting 
function [11] regular in all the cut w-plane and becoming infinite at m= +7. 


' See Marsnak [6]. 
* R{n} =the real part of 7. 
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On the other hand, for real and finite values of 7>1 which will be seen later 
to be the most important in the evaluation of (1), all functions h» are real, 
hn>0O, and also hy<1 at least for the first few values of m. One may there- 
fore guess, and it is also easily shown [11] that fy then has a purely imaginary 
pole w in the interval 0, +7, which is in fact the pole called above the main 
pole Wo. 

It is always possible to perform one of the integrations in (1) by the method 
of residues. Following WALLER [9], we integrate @ first, using this method and 
are then left with an integration in 7. The reverse order of integration can 
also be used [6], [11]. It is then convenient to introduce a variable 6 in- 
stead of 7, 


(8) 6=7—1. 
Then, for the main pole wp=@ (9), we put 
(9) cop (0) =i B (8). 


For real positive 6-values B(6) is always real and positive with B(0)=0. From 
the preceding it is also obvious that 


(10) lim B(6)=1. 


600 


A more detailed proof of this can be obtained by a direct study of B(@) for 
large 0 [11]. 

Concerning the other poles in the upper halfplane for real positive 6-values, 
they lie in general symmetrically around the imaginary axis with considerably 
larger imaginary parts than w,(6). For @+0co they all move towards the cuts. 
In performing the integration in w by closing the contour of integration by an 
infinite halfcircle in the upper half plane broken off by a path along the cut 
around the point +7 and forming the residues at the enclosed poles we see 
from (1) that the “secondary” poles will have very little importance compared 
to the main pole. In fact, we are going to evaluate the integral (1) only for 
such values of R and «x that the secondary poles can be neglected. 

We shall now consider more closely the main pole. In a small circle around 
the point 6=0 B(0) may be determined by a series 


(11) [B (0)? = 0? 0 (1 b, 0—b, 0 )=C?6(1 p10 Bo 0 Saye 


Here the constants C?, by, by ..., By, by... can be successively determined. 
They will depend on the values of h,, (7) and its derivatives at the point 7=1 
with the notation 


: 1! 
12 Hie alt 
( ) {n ad yf ih ()} 


taken for 7=1, (= symbol of differentiation, n=0, 1, 2,...). The expressions 
of some of the first constants in (11) in terms of h' have been calculated be- 


fore, obtaining [4]: C?=3hoh, ete. 
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However, we must have accurate knowledge of B(6) for much larger real 6 
than can be provided by (11). It can be obtained by analytic continuation of 
this expression of B(0). In fact, taking, in the cut w-plane and the 6-plane 
where (7) holds, those points @, 6 which make |f>1(@, 0)| limited to some 
positive value />!'(@, 6) is holomorphic both in w and 6 within these domains. 
Then, by the theorem of the permanence of the functional equation, the ana- 
lytic continuation w(6) of the expression (11) of @ (6) will always satisfy the 
equation f,'[w(6), 0]=0. In particular, the analytic continuation, along the 
real positive 6-axis and then along the central part of the path of integration, 
chosen sufficiently far from origin, always gives for w(6) the main pole «po (9) 
as no branch points occur on these paths. This makes our definition of B(6) 
more complete and also gives us a method of calculating it, facilitated for 
large 6 by (10). 

Calculating now the residue of fp at the main pole we obtain the value 


“3h 
(13) Res {fo (@, hoot = 32, G4) 
This defines a new function A(6) for which we have A(0)=1, and in a small 
circle around this point 


(14) A (0) =1+a, 0+ a 67+ a 07+ --- 


Just as B(6), A(0) may be calculated outside this circle as the analytic con- 
tinuation of the series (14). A direct study of A (0) for large real 6 shows that 


(15) lim A (0)=0 


6>00 


which easily follows from a similar relation derived by Wick [11]. The constants 
a may be successively calculated in terms of h’. However, once the 0, are 


calculated « may be obtained directly with the use of a simple formula 


0 
y=— aan (ho by RS Ph ecosiah: 
(16) a = — Ba Uh Be 0 ) 
For the derivation of this formula see Appendix A. It was given in an earlier 
paper [4] but without proof; for explicit expressions of the three first «, see 
this paper. 
Finally we obtain as the result of the integration in w 


16 


(17) F(R, BE) = 


Here the path of integration Lo is a straight line in the right half plane pa- 
rallel to the imaginary axis. (17) is obviously not exact, but the validity of 
the approximations made will be checked later. 
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3. Determination of the functions B(6) and A (@) 


The method we intend to use for the calculation of B(0) and A (0) was 


sketched above; for small |@| by means of (11) and (14) and for other 6 needed 


by analytic continuation of these series along the real positive 0-axis. 

In order to carry out this program in practice we must start with a thorough 
knowledge of B(@) and A(0) in the neighbourhood of origin, which can always 
be obtained by calculating sufficiently many coefficients in the developments 
(11) and (14). However, the derivation of explicit formulas of b, and « in 
h' and then substitution of the values of h’ is an extremely laborious procedure 
when » is not small. Therefore we have determined two recurrent systems of 
equations for the calculation of 6, and «, which were found to make the nu- 
merical work rather easy. /, is then calculated by means of (11). 

We define a function A,—1(0) by the continued fraction contained in (2), after 
substituting for w the value at the main pole 


2 


(18) An-1(O)= lini (OY + aoe [or (9) 22 “=)235m ). 
hn (O)] + 


An-1(9) is a regular function at least for sufficiently small values of |6| for 
which we write, in expanding hn_1(6), hn (0)... and inserting (11) 


(19) AO eee. (Gno=1) (n=1, 2, .. 
v=0 

The recurrent system is obtained from the relation 

(20) [An-1 (62) —hn=1(6)] An (9) = —7—p —[B (8) (w=1, 2, ...) 


written as an expansion for small |6|, observing that Ay (#)=0. The resulting 
system of equations for the a», and b, is best written as two formulas 


g “ ne 

¢ ~ haan n 5 

(21) hy Win » (cn aay i h — An, —p = 4 2 i CO? nee 
weal Le Hn —] aad 

me (n=2, 3,..., v=1, 2,9 

Wi 
sty olvesy a7) 
b= —1, C*=3hoh,, 
and 
4A hi 

26 = = 0 P 

(22) Oy—4 iad eee >: =A Q1,»—y (v= ils 2, oe eg 


/ reer a 
et la! ho 


Using (16) we then obtain by differentiation of (21) and (22) another recurrent 
system for the determination of the «,. For abbreviation we put 
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7) 


23 = hoy Ne: 
( ) Any = ho A (ha) (Any), 


and then we get 


y 


; ee 
(24) hy hn=1 by [ona An,» - a (Ana pare laa ) a, | ad 


u=1 ue! Rn—1 
—y" 
= ae | C Ay—1 (n= 2, 3 5 v= I, 2) 
With a=1, drno=0 (n=1, 2,...), 
v he 
(25) Cori) Whee be ae Cin (at ee te 
nae 


If we now want to determine all 6, and « up to a certain index »=k, we 
start with a calculation of all h' for which n+i<k+1 and then compute the 
hi hi 


corresponding and (n=1) and put them into a triangular scheme. 


segdinn 
a! hn 


Using (21) we then calculate an,1— 


A 
i! ho , 
7, ( 
nN 
to the left of a second triangular scheme, and then a»,1 and, remembering (22), 
b,. After that we proceed to the second column of the scheme, containing 


n=k), and put them in the first column 


Bn, 2 (n=k—1) and also get dy,2 and b,. We continue in this way until 


ji ha 
ha: 
the value of 6; has been obtained. The intermediate values of an,, are for 
every n value put in a row, beginning from the right. This makes the pro- 
hn-1 
v! Road 
consists according to (21) merely in applying the n:th row to the n—1 row of 
the scheme, forming the products with reversed signs of the contigous elements 
and adding the right member of (21). The calculation of such sums of pro- 
ducts is very rapid using an electric computing machine.’ The computation of 
the « is then made in the same way, beginning with all dn1 and a, dy,2 and 
% etc. and, in order to facilitate the formation of the sums of products in 
(24), in putting d,,, in a triangular scheme similar to that used before, as well 
as for each n value in a row just as was made with the ay,, above. 

The numerical calculations sketched above must be preceded by a computa- 
tion of the constants h’. Although exact explicit formulas may be derived for 


hi, remembering (4) and cy(y) in the form given by WatLrerR [9] this point 


requires further attention. For M-=1 this formula is certainly simple and h’ 
may be rapidly calculated in this way, but in the case M>1 a numerial com- 
putation is very involved, except for small values of m and 7. The difficulties 
arise mainly from the fact that the exact expressions are small with increasing 


cedure easy in practice; in fact the calculation of a certain value a@p—1,,— 


1 The machines actually used were Facit, mod. ESA, Atvidaberg, Sweden. 
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values of m and 7 but are given as algebraic sums of large terms, giving no | 
account of the right order of magnitude of h' for which, however, very simple — 


relations exist. For instance for carbon, if <24 or n=2 it is already difficult | 
to get accurate values of h'. Therefore, for M>1 we have derived simple | 


formulas suitable for numerical calculation (Appendix B). From these formulas 
it appears that for M>1 we need not take into account more than the first 
few ¢n(7) in (2), putting all others c,=0. This fact has been pointed out be-_ 
fore [6], but here we also may easily study in detail the validity of this ap- | 
proximation. Thus the numerical calculations show that for carbon c,(7) and 
for deuterium cg(y) are the highest coefficients cn(y) that we need take into 
account in (2). This approximation, in fact, is excellent for the calculation of | 
B(6) and A (0) even for very large real values of 6. It also follows that b;, 6; and — 


a; are all of the order of magnitude aa 


;. Thus the radius of convergence 09 


of the power series (11) and (14) must be roughly proportional to M. From 
the numerical results reproduced later in this paper we obtain as a very rough 
9° approximately valid for all values M and for 
the power series of both B(6) and A(6). In practice, using ten terms for ex- 
ample the series may be used for accurate calculation of B(#) and A (6) for 
values of |@| approximately up to 0.4 M. 

Now the analytic continuation of the developments (11) and (14) is performed 
with the classical method using Taylor series. The method is slightly modified 
in taking 1/0 as variable instead of @ and putting 


(26) y= 5 HW-B(), Lw)=4(’). 


Then the method was found to be very useful. 

The expressions of H(y), H?(y) and L(y) for large y-values then follow from (11) 
and (14). Taking now a real positive value yp so that 6) =1/y9<ge (something 
like 0.309), we calculate the values of H(y)* and L(y) and their derivatives at 
y=Yo and form the Taylor developments at this point. Taking a new point 
y, on the real positive y-axis, 0,;=1/y, > 9 but still a little smaller than go, 
the Taylor series at y, are calculated using those at yy. We now write the 
new series aS power series in y—y, 


estimation the value for 09 = 


(27) HT (y) = ho + hy (YY) + ha (y—m1)? +, 
(28) L(y) =Io +h, (y—y,) + Ia (y—ys)* + + - 
The constants ho, hy,..., lo, 1, ... are calculated in the way sketched above. 


The choice of only one intermediate point yp for calculating (27) and (28) and 
the rough positions indicated for yg and y, are facts mainly born out by the 
numerical calculations which will be reproduced later. In this way, (27) and 
(28) imply in any case an essential extension of the knowledge of B(6) and 
A (9). In practice, of course, only a limited number of coefficients in (27) and 
(28) may be calculated, and this number mainly depends on the number of 


1 Tt is somewhat simpler to start this calculation with the function H*(y). 
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b,, P, and a, which are initially computed. This also determines the accuracy 
of the values obtained for B(#) and A(@) for large real values of 0. This 
accuracy can be tested by putting y=0 in (27) and (28) because according to 
(10) and (15) the values H (0)=1, L(0)=0 are known. However, although (27) 
and (28) seem to show good convergence even for y=0, for very small |y| 
more accurate values can be obtained by taking developments around y=0, 
(@=co), and using (10) and (15). We now perform an analytic continuation 
of H(y) and L(y) to y=0, calculating H (0), H’(0)... L(0), L' (0)... using 
(27) and (28). But before we do that (27) and (28) are completed, if needed, 
with terms above those actually calculated, these rest terms estimated so that 
we get H(0)=1, L(0)=0 exactly. This procedure is of course a little arbitrary, 
but in practice these correction terms are small, applied only in order to in- 
crease somewhat the accuracy of the values obtained for the derivatives at 
y=0. We finally obtain in the neighbourhood of y=0 (@=©co) 


(29) H (y)=1 ky ky y” ks y? ‘Tis 


(30) Liy)=myy+mey?+---, 


where the new constants are calculated from (27), (28) as sketched above. 

Now we have explicit expressions of B(@#) and A(@) for all real positive 
§-values consisting for both, of three developments: around 6=0, at §=00 and 
at an intermediate point 9=6,. The question of the radius of convergence of 
the new series is best answered in detail in each special case, when numerical 
calculations have been performed. However, it seems in any case that the series 
at §=0 and §=co have a common part of convergence in the 6-plane. As 
we also have taken an intermediate point 0, then for any 6-value in this interval 
we need only take the first few terms of a rapidly convergent development to 
get accurate values of B(O) and A (6). 

As appears from (17) accuracy of the value of B(@) is more important than 
of A(6), especially for large real 6, as will be seen in detail later. Therefore most 
of the numerical work in this case has been devoted to B(6). In particular, it 
is convenient for the evaluation of the integral (17) to have some simple ap- 
proximate formulas of B(6), called B;(0), as well as the exact expressions. For 
small |9| we may for that purpose use (11), taking B,(0) as 


(31) B, (0) =C V6 (1—f, 8). 


For small |@| this is a good approximation, and for any real positive 6-value 
when this is not valid, it was found numerically that another simple expres- 
sion could be used 

b 


2 ,(6)=a— ——. 
(32) B,(9)=a——" 


Here, a, b, and ¢ are constants to be determined numerically. However, if we 
want great accuracy, we cannot use the same values of the constants in the 
whole interval from a small 6-value to @=0co. But in dividing this interval 
into only a few smaller intervals, and taking, for each of these, certain values 
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of a, b,c, (32) was found to give a surprisingly good approximation not even 
of B(0) but also, when differentiated, a good approximation of the first. deri- 


vatives of B(0). In particular, according to (29) for sufficiently large real 0 
we may take a=1, b=k,(>0) and c= —hy/ky. For the remaining intervals a _ 


is in all cases close to 1, 6 and especially ¢ vary more (b>0). 
Using the reverse order of integration in (1) Wick [11] has outlined a pro- 
cedure for calculating the functions then corresponding to B(@) and A (6) for 


real positive 0. For small 6-values he uses developments, too, but taking only 


some few terms of them. For intermediate 6-values he sketches a purely nu- 
merical method, giving the values of the functions only point for point instead _ 


of as an explicit expression. Finally for large 0-values Wick attacks the problem _ 


by examining directly the functions in the neighbourhood of 6=00, and de- 


scribes a method of calculating k,, k, and m,. However, the direct study of 
the expansions (29) and (30) around 6=00, is difficult because no analytic 
formulas of the coefficients can be derived and the computation of ky, for ex- — 


ample is involved already. §=co is therefore not a convenient starting-point 
for a close investigation of B(@) and A(6), in sharp contrast with the point 
6=0, from which in fact we have calculated these functions for all 6-values 


in this interval. The separate methods developped by Wick respectively for — 


intermediate and large 6-values are not necessary. Besides, the method sketched 
here gives analytic and, except possibly for large 0, more accurate expressions 
of B(6) and A (6). 


4. Evaluation of formulas for the neutron density 


It is now possible to perform the integration in 0 which was left in (17). We 
shall use the saddle-point method, applied to this problem before by WALLER [9] 
and by Wick [11]. The integrand in (17) consists of the product of the slowly 
variable function A(0#) and an exponential function, being rapidly variable com- 
pared to A(O) provided that not both w and R are small. We may deform 
the path of integration LZ in the 6-plane to pass through a saddle-point along 


a path of “steepest descent”. Along this path in both directions the real part | 


of the exponent is rapidly decreasing from a largest value at the saddle-point, 


the neighbourhood of which is the only part of ZL that gives an essential con- | 
tribution to the integral. This may then easily be evaluated by developing | 


the integrand around this point [2]. In our case we have only one saddle- 


point lying always on the real positive 0-axis. Now we want to get the integral | 


as an analytic expression in # and R and so, to begin with we must be able 


to express the saddle-point analytically in 2 and R. The method of integra- | 


tion, as it has been sketched above, is therefore slightly modified in the 
following way. Instead of the exponential function in (17) we take as the 
rapidly varying function that with B(6) substituted by B; (0), and the difference 
B (0)—B;(0), being small may even when multiplied with —R in the exponent 
be considered as a slowly variable function just as A(0). We thus determine 
the saddle-point 0; by the equation 


(33) a—RB;(0)=0, with the root 0= 6s. 
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Thus 6; is a function of x and R only in the combination R/a, a fact which 
has been stressed before [11]. B’(0) and By (0) are real, positive, monotonically 
decreasing functions when @ varies along the real positive 0-axis with the limit- 
ing values 0, for 6-co. Then 6, must be real, positive and monotonically in- 
creasing with R/x with lim 0,=0o, for R/x>0o, (6;=0, for R/x=0); this is 
seen in detail from (31) and (32). The division of the real positive 0-axis into 
intervals with different values of a, b, and ¢ which was found to be necessary 
in order to get good approximations B;(@) then simply involves a similar divi- 
sion of the values of R/a into some few intervals. 

For the evaluation of the integral (17) we should now develop the integrand 
around §=6; in the usual manner and then integrate, the result being a cer- 
tain series development for F)(R, #). However, we must first observe that 
B(#) is not regular for 6=0, thus for R/x-0 ie. 6;>0 this method cannot 
be used. Therefore, at least for small values of R/x, (17) must first be trans- 
formed by substituting 


(34) 6=q? (+ VO=9Q). 
Introducing then for simplicity 
(35) S(y)=B(¢"), U(p)=~A(¢), 


(17) is transformed into 


= Q 3hy I { x2p—RS(p) 
(36) CNC HP aa dg U (y) e* uy 
L 


P 


where the path of integration L, originally is the picture in the g-plane of the 
path of integration J» in the 6-plane in (17). Now (36) may be evaluated in 
quite the same way as is outlined above for (17) and we observe that (34) 1s 
valid for the saddle-point ys in the g-plane, i.e. gs= +V6,. The explicit values 
of 6; and gs are 


(37) Gs = ———,, with t= 


ETS ni Pe 
_ Re pegs 
e(la az) ‘ Z 


for small (9s) within the range of validity of the approximation (31) Le. 
small R/x, and 


(38) angi-|/* 


for larger values of R/x and 0, where (32) is valid. 0;=gs<1/6, may be taken 
20 aie 
as an approximate upper limit for the use of (37), giving R/a#< G / os 

1 
for larger values of R/ax (38) is the best approximation. (37) was derived and 
first used by WALLER [9], it has also been used by the present author before 
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[3], [4]. In the special case R/x—>0co, when b=k, and c= —k,/k, (38) has been — 
used earlier by Wick [11]. | 

We now form the Taylor developments of «y*—RS(p) and of U(q) at | 
~=9s, taking ! 


(39) p= Pstit. 


All terms of the exponent, except x y;—RS (ps) and —[a—RS" (gs)/2]¢? are then 
expanded and we integrate from t=—co tot=-+oco. The path of integration is 
thus certainly not the correct path of “‘steepest descent” but only tangent to this — 
at the saddle-point; however, even along the path chosen here, the real part 
of the exponent decreases sufficiently rapidly from the point g=qgs and only 
a small central part is of importance. In this way we get from (36) an ex- | 
pression for F(R, Z), having the form of an infinite double series containing © 
x, R and the functions S and U and its derivatives at p=, and being valid — 
for any value of R/az. But when (38) holds we might as well use the original — 
integral (17) for the calculation of Fy (R, #), in fact, the introduction of the 
functions S and U and the use of (36) is then an unnecessary complication, 
the derivatives of B(#) and A(@) being easier to calculate than those of S(q) 
and U(qg) in this case. (17) is evaluated in the same way as (36) above, taking 
even here as path of integration a straight line through 0=6, parallel to the 
imaginary axis. The expansion to be integrated is of course of identically the 
same type as that above in (36) and we may therefore write both expressions 
for Fy (R, £) in one formula 


Q-3h, i | 15ni 3n,. 3465 ni 
40) F(R, E) = ——— @%%— RBs) k. M, 
ne poli) (2Vn)®RE~ iia Upc ia 512 


_ 945 ni ne 7 105 m3 L05nyng 15 mq u 3n, 31573 
na! hee ey) ih oe a ee 64 


105m n. 155 Li, L056 j.. HLS tas 
16 8 + loa | 32) 8 soak 


+ tg | SP + | a; mli+-| peal. 


Here, for abbreviation, we have 


u j 
(41) M; = >, mii (é=0, 1, 2,...). 


j=0 


In (40) and (41) we have introduced hk, e, m; and n; as notations for certain 
functions of R, « (and M) (e is in general a small quantity). For all values 
of R/x we may use 
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ley aloe Pg Co] Ne eee ore 
2) =2)0— SS" (p)] » €= 5[2rps—RS' (ps)], 


~_ 


k\i VU (qs) ; k\i RS® (, 
m= (5) peda lien ie Reree UY m2 = (5) _f (¢=3, 4,...), 


obtained from (36). However, as was mentioned before, (42) is convenient to 
use only for small values of R/x when (37) is valid, otherwise it is simpler 
to take 


bf 


=" 
(43) Ee |- ty (6 | ERLE (8) 
(i) (i) 
: ne (7 =), 1; 2, oe Ns ri ag ee eae (=3, 4, Ar a) 


mz = ki — : 
a! a! 


these values being obtained from (17) and used when (38) holds. Regarding 
M; and n; as small quantities of the 7:th order of magnitude, in the expansion 
of F)(R, E) in (40) terms up to the fourth order have been calculated, and 
the formula therefore looks rather complicated. It is well known that, in general, 
a series of this type is not convergent but still can be used for practical cal- 
culation if the first terms decrease sufficiently rapidly. We may also terminate 
the series earlier if the first terms in (40) are small, this of course also depends 
on the accuracy wanted. 

The formula for Fy)(R, #) given here contains a great many earlier formulas 
as special cases. In neglecting all terms but the first in the expansions (11), 
(14) and (40) we obtain the age-approximation of the neutron density. The con- 
dition for R/z in this case is, according to (37) R/x<1/CV3£,. For larger 
values of R/x, approximately up to this limit the formula of WaLier [9] can 
be used. This was obtained by taking one more term in the series (11) and 
(14) and using (37) and is contained in (40) by taking terms up to the second 
order. For still larger values of R/x more terms must be taken into account 
in (11) and (14) and the formula becomes more involved. For all R/ax ap- 
proximately up to 2V2/CV3,, ie. as long as (37) is valid, an earlier formula 
by the author [4] may be used which also contains the development of Mar- 
sHAK [6]. In the derivation of this formula only the main part, (#t)~“?, of : 
was kept in the term —4/?/k? of the exponent before the integration and in 
the expansion the explicit dependence on R, x and the coefficients of the series 
(11) and (14) was written out, taking these into account up to /3z and a3. It 
may be pointed out, however, that the way of writing the expansion of Fy (R, £) 
in that paper is not strictly consistent for the largest R/a in question when 
(40) is somewhat better. This is on the other hand not a strictly consistent 
development for small R/xz, but it seems unnecessary to rewrite it for this case, 
when some of the simpler formulas mentioned above may be used. Finally, 


the condition for these formulas and for the expansion (40) to be useful (when 
R/x = 2V2/CV3£,) may be approximately written at>1 (say 210), obtained 
from (40), (42). 


Or 
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When R/z surpasses the above mentioned value, (38) gives the best value 
of the saddle-point and (43) is simplest for the terms of the expansion of 


F\(R, E). Then, for intermediate values of R/a (40) now contains just as be- ; 


fore the numerical formula by Wick [11] and for R/xz+oco his asymptotic for- 
mula, obtained by taking only the term My of the expansion (40). If the ap- 
proximation B,(0) of B (0) (32) is determined sufficiently well, we may even here, 
instead .of the accurate but rather involved formula (40), derive a simpler ap- 
proximate formula, obtained using B;(9) instead of B(6) in (40) and (43). We 
then get 


: 
16Vb Rx 


Q -3 hy ba 
(2 Vx)? EB (Ra) 


— set | ood F 15 + |- 
512bRx : e 14 64VbRez 


A" (0s) a = .- | -AeOV2r [+ 
ey \ ai4  @64VbRe x 2 |16— 


ET Bee [peed 


x J 


(44) F(R, £) = e~tR+2V bRz—ex 4 (9s) [: 


(iv : 
ae is decreasing for increasing 7, R/x and 6; which may be seen from 
a! 
(30). Then it also appears that (44) and (40) are essentially expansions in powers 
of 1/VbRaxw. The condition for these expansions to be useful, for the values 


of R/x in question, is therefore VbRa>1, (say 24). In (44) a, 6 and c are 
constant within the intervals of R/z mentioned above. The fact that b has to 
be altered at the endpoints of these intervals essentially determines the ac- 
curacy of (44) compared to (40) for intermediate R/a and gives then a maxi- 
mum error of (44) of a few %. For very large R/x the accuracy of both (40) and 
(44) is rather critically dependent on the accuracy of the values of B(6;) and 
Bs; (9s) respectively, being ultimately restricted by the limited knowledge of the 
values of M. In this case we must, therefore, always be content with a good 
value of log [F9(R, #)]. This is realised by (40) and (44) even though B(@) has 
not been so accurately determined in this case which seems to be of less import- 
ance. 

Finally we must discuss the validity of (40) with respect to the approxima- 
tions made in deriving (17) from (1). In the first place the contribution to Fo 
in (1) from all other poles of fp than the main pole was neglected. The reason 
for this approximation to be legitimate was touched a little above; more clo- 
sely it has been investigated by Wick[11] and his results were also confirmed 
by rough numerical calculations which were made for the ‘‘h,-pole’”. Secondly, 
the integral around the cut +7 —-+ 700 in the w-plane was neglected. This 
“branch-point” contribution is also of a type studied before both in diffusion 


theory [7] and slowing down theory [6]. Therefore it is not necessary to make | 


a detailed calculation here. The essential point is that this contribution is at _ 


most the magnitude e-” and is therefore neglible (cf. (44)) except for very small 
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R when also w is small. Summing up, the conditions given above for the ex- 
pansions of Fy) (Rk, £) to be useful, are also sufficient to ensure the validity of 
the approximations made in deriving (40) from (1). 


5. Applications to special values of M 


Now we are going to present some numerical results for different M values 
and also study in general how B(0) and A (0) depend on M. Detailed numerical 
calculations have been performed for deuterium and carbon. Furthermore, for 
the coefficients of (11) and (14) which are the most important of the three 
developments of B(#) and A(6@) and the only ones having analytical expressions 
in M, explicit formulas, for M@+}1 are given for some of the first coefficients 
in form of power series in 1/M. For completeness and comparison the corre- 
sponding values for M=1 were also calculated, even though this has little 
practical interest on account of the assumption of a constant mean free path. 
We first present these general values for M>}1, followed in parentheses by those 
for M=1. The constant h,, occurring for instance in (40), has the value 


h, =1 
eat FAT 
foes ) 
ome V/ ae 3M -+ oy | for M>}1, (C=1 for M=1), 
Sb 1 Che alt | (30) 
a M 3M 45M? 30)” 
2.870 0.844 0.59 
(45) 3.057 0.615 0.65 
fa = =a (  M +o} eee eal 
1.82 0.24 
i vw ( va -) (1.238), 
0.72 
ps= we (2-01) 
3.867 0.782 0.43 ( 4 
eee, eT | Sell Ahatli =) 
: M M M? 5 
(46) | 
8.008 0.863 0.50 
he — We (1 M + Me -) (0.3606), 


1 For more complete numerical results for this case cf. Wick [11]. 
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419 0.43 
ma 207 (1 LE ule ) (0.3520), 
4.95 0.55 ss 
(46) | %4= yea (: iw) ‘hae 
4.7 
ee aie (3.01). 


The developments with the three terms calculated in 1/M give good values 
even for deuterium, and for carbon the agreement is excellent; the remaining 
series give in any case even for deuterium at least the right order of magni- 
tude. For a detailed study of the accuracy of (45) and (46) we refer to the 
direct calculations for deuterium and carbon which will now be quoted. 

The coefficients calculated in (45) and (46) are sufficient for the knowledge 
of the neutron density within a wide range around the source and in any case 
for all values of R/a when (37) is used. However, for an accurate performance 
of the analytic continuation more coefficients are needed. It would be quite 
possible to derive general formulas in the form of power series in 1/M even 
for higher coefficients but it is rather laborious. Besides, in performing the 
analytic continuation of (11) and (14) and using (26) the analytic dependence 
of M will in any case be lost so that the coefficients of (27), (28), (29) and 
(30) must be determined separately for each M-value. Therefore we refrained 
from computing more f, and a, values in the general case, restricting us to 
deuterium and carbon. 

Deuterium. We have taken M1=1.9963, yo=3, 09=1/3 and y,=6,=1. The 
values then obtained for the coefficients of the developments of B(@) and A (6): 
(11), (27), (29) and (14), (28), (30) are given in Table 1 taken in this order 
with the values referring to B(0) im the three first columns. 


C= 1.2046 pobre 
mH B n h n Ky, | tn | ly My 
0 (—1) 0.7785 (—1) | | 0.399 (0) 
1 0.63273 —(0.1540 0.337 = ee 0.255 0.61 
2 472350 0.0432 —0.21 | 12769) 00087 —0.3 
3 0.29596 | —0.0137 0.2 | —0.9254 | 0.033 
4 —0.138654 | 0.0052 0.4430 | --0.010 
5 0.03379 — (0.0025 —0).0793 0,003 
Ga | 0.00647 —0.0531 
7 —0.00700 | 0.0166 
8 —(.00333 0.0502 
9 0.00745 | — 0.0569 
10 —(0.00361 0.0090 
11 —0,00185 | 0.0348 | 
12 0.00364 | 
13 =0.00155 | | 
14 —0,0012 | 
15 0.0020 
16 — 0.0007 | 


* The values of M are obtained from Berar, “Elementary Nuclear Theory”’. 
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In this table the last figure of the values of kn, my, and also of the four 
highest hk» and /, has to be regarded as uncertain since it has been obtained 
by estimation of the rest term of a series. We have, however, possibilities of 
an indirect test of the values for instance by comparison of the series within 
a part of the 6-plane of common convergence, or for (27) and (29) by putting 
y=0. These tests give positive results. Concerning the long columns of coef- 
ficients P, and a, of course only the first are needed explicitly in (40), at most 
five or six, the others being calculated only in order to obtain the values in 
the other columns. As appears from these, the new series are rapidly con- 
vergent and give good values of H(y) and L(y), (27) and (28) for |y—1]|< 0.8, 
(29) and (30) for |y|<0.2. Thus in (40) it is best to use for B(0,) and A (9,), 
(11) and (14) for 0,;=5/9, then (27) and (28) up to 0; <5 and in the remaining 
interval up to co (29) and (30). These values of 0; correspond to the following 
values of R/x: 2.97, 92 (and oo), obtained from the values of a, b, and c in 
(32) and (38) which will now be reproduced in Table 2. This table also con-. 
tains the values of certain constants d, e and f introduced in an approximate 
expression of A (6), As;(9), similar to (32) 


(2 


(47) A, (9) = d+ 5 


This was found to be a good approximation and in (44) we may use (47) for 
A(6), thus getting a somewhat simpler formula. In Table 2 the two first co- 
lumns give the approximate intervals of 6; and R/z corresponding to unchanged 
values of a...f so that each row of values a... f is put between the ap- 
proximate endpoints of the intervals within which these values ought to be used. 


Table 2 
0, R/x a b G d e f 
4 

ee me 0.9871 0.2845 | 0.3637 | —0.047 | 0.762 0,712 

70 nae | 0.9971 0.3126 0.430 0 0.612 0.53 
10 | 330 1 0.3325 0.534 0 0.612 0.53 

| a 1 0.337 0.62 ) 0.61 0.5 
co CO 


The values in the last row were obtained directly from /, and m, in Table 1 
with b=k,, c=—kp/k,, e=m, and f=—me,/m,, the three first of which are 
contained in the asymptotic formula for R-~co of the neutron density. The 
other values were calculated simply by taking B(6) and A (0) respectively at 
three points in the corresponding intervals and give a maximum error at any 
point of about 10* and 10° respectively. 


Carbon. We have M=11.905, yp =0.6, 69=5/3 and y,=0.3, 0,=10/3. The 
values of the coefficients are given in Table 3. 
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Table 3 
OC =0.67094 
eS Sl i ea a 
n B,,° 10" hy ky a,,* 10% Ly My 
0 (—1) 0.8326 (—1) (1) 0.406 (0) 
1 1.4303 — 0.5406 0.49 — 3.0435 1.03 183 
2 — 1.8890 0.2997 0.4 5.2583 — 1.44 2 
3 1.7255 0,324 —(,4 — 5.2038 a 
4 — 0.8984 — 1.59 2.4295 2.0 
5 0.2097 3.2 — 1.0422 ho 
6 = (VSG — 23 3.7191 
7 0.7172 =i — 3.324 
8 - 0.1526 6-107 — 6.413 
9 2781 == Doi 8.301 
10 = 0.2197 6-10? 16.60 
ll 2.927 — 25.6 
12 = 0029 
13 — 6.895 
14 1.407 
15 16.62 
16 = 6:89 
17 —40.49 
18 26.10 
19 98.5 
20 — 89 


What was said about the values of hy, kn, In and my, for deuterium in Table 1 
is also valid here, and in particular the values of k,, kz and mz are very rough. 
The developments (27) and (28) show poorer convergence at y=0 in this case, 
than for deuterium and therefore the analytic continuation of these series to 
y=0 was combined with graphical interpolation for the calculation of k, and 
mn. Rapidly convergent expressions of B(0s;) and A (Os) in (40) are obtained 
if we take (11) and (14) for 6;22.5, then (27) and (28) up to 6;59 and in 
the remaining interval up to co (29) and (30). The corresponding values of 
R/«x are 13.1, 111 (and co), obtained from Table 4 which contains the values 


of a...f, arranged in the same way as Table 2. 
Table 4 
SS - _ — 
O, | Ria a b c d e f 
1 ; 
‘ ae 1.0683 1.0574 1.1947 | 
ee a 1.0421 (0.8938 0.9338 | -0.187 3.761 3.01 
6 | a 1.0089 0.6164 0.148 —0.055 2.203 1.43 
13 300 1.0029 0.548 —0.19 | —0.004 1.48 0 
Sire eee 1 0.49 5 colhlied 1.3 —2 


These values were calculated in the same way and give about the same ac- 
curacy for B(#) and A(@) as those for deuterium in Table 2. Finally as an 
illustration the logarithmic neutron density corresponding to E)=0.1 Mev, 
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E=1 ev is shown graphically in Figure 1 together with curves obtained from 
the age-approximation and the asymptotic formula. 

We have now collected enough information of the functions B(6) and A (6) 
to be able to study in general their dependence on M, at least if only the 
values of physical interest M=1, 2,... (approximately) and real positive 6- 
values are considered. If we write for a moment By (0) and Ay (0) (45) and 
(46) imply at least for small constant 6-values By, (0)>Bm,(0) and Aw, (0)< 
Ay,(9) for 25M,<M,. These relations also seem to hold in the whole inter- 
val for 0<6@<co but then only for sufficiently large values of M, and Mg. 
For large 6 the relations are valid because then m, and k, increase with M; 
in fact, in the limiting case Moco for which some preliminary computations 
were made, taking 6/M as variable we get k, and m, proportional to M which 
justitifies our statement, at least for large WM. But they seem to be valid for 
a value of M as low as 5 (roughly), though we would need more numerical 
calculations to study this im detail. However, for the smallest M-values we 
have on the contrary B, (9)<B, (0) and A, (6)>A,(6) for constant 6, 0<@<co, 
Thus. of all physically possible functions By (9) and Ay () those for deuterium 
are the largest and the smallest respectively for small 6-values; for larger 6- 
values the functions in an interval 2<M/=<6 (roughly) seem to lie close to each 
other and the extreme values of By and Ay seem to be attained within this 
interval.? 

More detailed information on the dependence on M for large 6-values is ob- 
tained from the calculations for deuterium and carbon, being characteristic re- 
presentatives of the case of a small and large M-values. It would perhaps be 
expected, that the analytic continuation would be easier for carbon than for 
deuterium as then the original developments (11) and (14) have much larger 
circles of convergence. But a comparison of Table 1 and Table 3 shows that 
for large # it is just the contrary, as we have pointed out above also. For 
example, we get poorer values of k, and m, for carbon than for deuterium. 
These unexpected facts are best clarified if H (y) and L(y) are studied graph- 
ically for small y. For deuterium these functions behave gently and reach their 
values at y=0 monotonically and with no change of the sign of the curvature 
for decreasing y. But for carbon there is a point of inflexion in the neigh- 
bourhood of y=0 which causes trouble when the analytic continuation is per- 
formed beyond it. This behaviour is also seen from Table 4 which shows some 
characteristic differences from Table 2. The values of a... f vary much more 
for carbon than for deuterium. In particular, for carbon ¢ (and /) gets nega- 
tive for large 9, which is the most direct indication of the existence of a point 
of inflexion. 

What has been said above for carbon is certainly also valid in general for 
large M-values. In fact, the preliminary computations in the limiting case 
M-co show that these properties of B(6) and A (9) are then still more marked. 
But this also has the effect that the asymptotic formula for R->co of the 
neutron density is a much poorer approximation for intermediate R/a for large 
M than for small M. It is therefore also of less practical importance in the 
former case to perform the analytic continuation accurately as far as to y=0, 
determining nice values of the constants of the asymptotic formula. 


’ Cf. also for ky some computations by Marsuax [6]. 
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6. Discussion of the formulas for the neutron density 


We have now derived general formulas for the neutron density and studied, 
in general, the dependence on M of the included quantities. Therefore we are 
finally ready for a discussion of the neutron density and its dependence on the 
three parameters R, x(#) and M. Then we write for (40) 


eh, 
(48) Fy(R, £) = Byatt? sa 


where v=vy(R/a) and also, to the first approximation, uw =u (R/xz). The most 
important of these functions is obviously y. It seems worth mentioning first 
that » very easily can be obtained graphically. In fact, let us select in the 
diagram of B(@), as a function of 6, out of all straight lines having a slope 
a/R, that line which is tangent to B(6). Then obviously the point of tangency 
has the abscissa 0; (if the small difference B(0;)—B,(0s) is neglected), and the 
tangent cuts the ordinate axis at a point where the ordinate is precisely equal 
to v. But we have also derived analytic expressions for the function vy =vy (R/2). 
When R/x varies from 0 to co y varies monotonically from 0 to 1. At first 
we have v=ROC?/4a, thus » varies rapidly with R/z and for large M is roughly 
proportional to 1/M (v+3R/2a”M). This holds, however, only within the small 
range of validity of the age-approximation or roughly when R/x<1/C V3 py 
(=M/6 for large M; ~0.6 for M=2), giving y<1/4. For larger values of 
R/x in an intermediate range, vy changes to a more slowly variable function of 
R/a which also becomes more involved. At first » transforms into the form 
derived by WALLER[9] v = eee are 

4¢ 3(1+7)* 
curacy with more terms when R/« increases [4]. Then, roughly when R/z surpasses 
the value 2V2/CV3 By (= M/2 for large M) and y attains the value 0.4 (roughly) 


; , : 64 ~<6% : 
we obtain again a simpler form of y=a—2 | —~ + —~ which becomes more 


(for t cf. (37)) completed for ac- 


and more slowly variable with increasing R/z. It is interesting to notice now, 
however, that the slowest variation with R/z soon occurs for small M in con- 
trast to the case when R/z is small. This is also connected with the fact born 
out by the detailed study of By (0) above, that the asymptotic form of y [11] 
(a=1, b=k,, c=—k,/k,) may be used as an approximation earlier, ie. for 
smaller R/x, for small M-values than for large M-values. In general, if we 
consider » for any given value of R/x (+0, +00) it follows let us say from 
the graphical method indicated above, and the discussion of By (0), that vu 
decreases when M increases and M is sufficiently large (in any case for M=6). 
For smaller values of M down to M=2 the conditions are different for small 
and large R/x but in any case the values of vy for these M lie close together. 
However, we always have »,>, for 0<R/a<oo, 

Though the term e~’” in (48) is the dominating one we may also for com- 
Pale Vzx)8, this term starts 
with the value xC (roughly proportional to 1/VM for large M) for R/x=0, 


pleteness discuss the term yw. Taking x = (s= 
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and is at first almost constant for very small R/z. But very soon outside the 


. ' (roughly) [9], [4] and now 


age-approximation it transforms into j= C—— 
Vr(l+t 
decreases more and more rapidly. When R/x has increased to the value 
2V2/CV3,, (roughly), we have =x" (a/ RY A (6s) (approximately), varying 
at first roughly as 2/R but for increasing R/x transforming into the asymptotic 
value w= xm, (ky)~“*(x/R)* [11]. This quantity is increasing with M for large 
M so that the dependence on WM is reversed compared to the value jy (0). 
But of course this fact is of little importance in (48) compared to the influence 
of the exponential. 

This discussion of » and w also gives the full picture of the neutron density 
as function of R, x and M with the restriction of course that the conditions 
must be fulfilled, for (40) to be valid, ie. roughly «= 10 for small R/ax and 


VbRx=4 for large R/x, these are also essential for the discussion of i above. 
We may study for instance, for every value of WM, the space distribution of 
neutrons of constant energy, i.e. constant x. However, it is only for x=10, 
approximately, that we get the complete space distribution. For higher neutron 
energies, 1.e. smaller xz, we have a small dominion around the source, extending 
with decreasing x where our formulas become uncertain but outside which we 
still have the complete space distribution. This may be described for inter- 
mediate R already as a quasi-exponential decay with a decay-length >1 slowly 
decreasing towards the limit 1 (=A). 

Let us then study the energy distribution of neutrons at a certain distance 
R. The dependence of x for the factors w and e-’” is then described by the 
discussion above. The factor left in (48) is decreasing with increasing # except 
when £ is close to Hy but otherwise all factors co-operate to make the energy 
distribution function decreasing with increasing #. This also makes the restric- 
tion, for the complete space distribution to slow neutrons pointed above, less 
important. When both z and R are small, source neutrons and neutrons having 
suffered only some few impacts are also of importance and our formulas are then 
insufficient. As an illustration of what this restriction implies for the energy 
distribution, let us take for carbon R=20. Then it seems possible to use the 
formulas roughly for z>1 which indeed is a very modest restriction. In general, 
at distances of this magnitude, still small enough to be of practical mterest, our 
formulas give almost the complete energy distribution. 

Finally, taking constant values of R and x, we may study how the neutron 
density depends on M. Taking first R=0 and constant values of Hy and FL, Fo 
increases rather much from M=1 to M=2, varies then slowly for the next 
M-values (maximum for M=4) and decreases then roughly as 1/V M. But for 
increasing R this picture is soon reversed with a slow variation and a minimum 
in the interval 2<M<6 and then Fy increases with M for larger M-values and 
also when we pass to M=1. However, in making this comparison between 
different moderators it is important to remember that in general we have dif- 
ferent length units for them, these being from the beginning chosen as the mean 
free paths Ay of the neutrons. Thus we have for each M R=r/Am, if r is 
the distance from the source measured in ordinary units. To make this com- 
parison complete we must therefore finally apply a spherical dilation of the 
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100 200 r (cm) 


Fig. 1. Log Fy + constant, calculated for carbon from: 


formula (40) (middle curve) 
age-approximation (lower curve) Ey=0.1 Mev, H=1 ev, 
eeempeate formula (upper curve, very rough) w#=11.51, A=2.:65 cm. 


R-space for each M, the scale of length being Ay to 1, and then we obtain 
Fo(r, EZ). A full discussion of the dependence on M for the neutron density 
should therefore also require a consideration of the variation of Ay with M 
which, however, is beyond the scope of this paper. But it is worth mentioning 
that this has a particularly large effect for hydrogen (in water for example) with 
its very small mean free path (EZ) small), making it a most effective element 
for purposes of slowing down, much better than say deuterium, in spite of the 
fact, pointed out above, that the slowing down per mean free path is better in 
deuterium. 

Finally concerning the practical applicability of our formulas it must be 
stressed that the assumption of constant mean free path for energies in the 
interval H—Ep, is justified only for sufficiently small #y. For carbon this holds 
for Hy <0.1 Mev', with 2=2.65 cm (taking the carbon density as 1.60 g/cm’). 
For deuterium [6] the approximation is good for Ey <= 0.7 Mev, the relative varia- 
tion of 2 then being less than 5% (roughly). An extension to the more im- 
portant case of heavy water will be given in a later paper. For hydrogen and 
water [6] the approximation is justified only when Ep is at most of the order 
0.01 Mev and the formulas therefore seem to be of very little practical use in 
this case. 

_ For higher values of Hy we may try to use in our formulas some mean value 
A of the actual values of A within the slowing down energy interval. Then of 
course the formulas can not be expected to give the same accuracy as when 
Hy 1s small. However, in chosing J judiciously we may still hope to get useful 
formulas, at least qualitatively, for moderate values of Ey. Considering A as a 


function. of « A=A(zx), age-theory [6] gives us the information that we should 
choose A as 


Ne 7 a 
Cf. [3] and [6], where also further references to measurements of A are given. 
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(49) a (é) dé, 


™I 

I 
So ir 
— 


for neutrons of energy E=E,e~*. When the conditions given above, for the 
age-approximation to be valid, are fulfilled, i.e. small R/ax and large x (49) must 
be a good value of 4 in our formulas. Supposing now the most frequently occurring 
case of a mean free path which decreases with decreasing #, varying between Amax 
and Amin, (49) gives a value of A which is not much larger than Amin (if not 
Amax/Amin is too large) and varying slowly with x. Now one must expect, for 
physical reasons, that the fast neutrons corresponding to A=Amax increase in 
mportance with increasing values of R (and EF); mathematically this has also 
been shown by Wick [11] in his formula of the asymptotic neutron density for 
this case. Therefore, (49) gives too small a value of 2 for increasing R (and 
decreasing x). In fact, for R->oo we must have yR=yr/A>r/admax' [11]. Thus, 
in order to make our formulas more useful, instead of taking a mean value /, 
constant for all r and x, we start with the value (49), variable with x, for 
r=(0 and then for each z-value we let A slowly increase with r towards the 
limit Amax so that »/A is monotonically increasing towards 1/Amax and rv/A> 
+r/dmax- With that prescription our formulas are useful for small R/z and also 
give the right main term of the logarithmic neutron density for r>co and 
therefore it does not seem too bold to hope that they will be at least qual- 
itatively applicable for all distances for a moderate variation of A. For inter- 
mediate R/x the slow increase with R/x of » in (48) was stressed above. 
For the coefficient of —r in the exponent v/A the slow increase of 2 means 
that this factor will be still more slowly increasing. This makes the quasi- 
exponential behaviour of the neutron density for a constant value of x (in 
disregarding the factor 1/r (roughly)) which was pointed out before, still 
more prominent for intermediate values of R/z; the decay-length is now larger 
than Amax and slowly decreasing towards dAmax (cf. also BorHE[1]). This be- 
haviour is also born out by experiments, by JENSEN and RITTER among 
others [5]. 

Though it seems possible in this way to reduce the error introduced by the 
assumption of a constant mean free path it is in general necessary for high 
values of Hy to take into account the energy dependence from the start if we 
want accurate formulas. Some important results for this case are due to 
Wick [11]. In a previous paper the present author derived formulas of the 
neutron density at small distances from the source, assuming a special form of 
energy dependence for A [3], cf. also [4]. In a later paper this problem will 
be taken up again and more general formulas given in order to complete the 
present paper. 

In conclusion we point out the slight modifications necessary for taking ab- 
sorption into account, assuming the scattering mean free path /s, as well as 
the total mean free path A, to be constant. It implies a multiplication of q 
and all ¢,(7) in (2), (4) with the constant *=A/As(0<%<1) [8], [9]. The main 
Phaenae © Uyare WOLO NO a, is the root of the equation BO) = 1 OR pee 
At the same time also / has to be corrected [11], Hotre (to be published in Ark. Fys.). 


™ More accurately: r// 
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pole @ is then connected with the root <1 of the equation ho (7)=1— 
—x0¢y(n)=0. mo may be determined numerically for a given value of x. Taking | 
0=n—no instead of (8), the functions w (9), B(#) and A (6) are defined and 
may be determined exactly as before, with the difference only that the constants 
h’ now refer to the values of the functions hy (y)=1—xen(y) and their deri- 
vatives at 7=1. Substituting only (Zo/Z)™ for the factor Z/# and multiplying 
with x, the neutron density will have the same form as (40) and (48), with 
the same general properties of the functions ~ andy. For any value of R/x#+0, 
y, in particular, will be increased compared to the case treated before (x=1), | 
and the neutron density, of course, is reduced for all R and zw. 


Appendix A 


Proof of (16). Introducing for simplicity a new variable 


(50) s=0, 59 (0) = 08 (6), 
we write 

1 : | S890) 
51 Res { \ bsp anaes Oey 
(51) es ilo (9, CO) 1 NS) 2 wo (8) Sree f-1(0, s) 


For small values of 0 and s we have 


Camm: ‘es | 
9 ho eae SA wea 


(52) 1-1 (0, jo 9 ae 


and, remembering the value of C?=3hoh,, from (9) and (11) 
(53) 89 (8) = —3 hy 8 (ho—ho by 9—ho by O?— - --). 


The functions (52) and (53) also depend of the constants h'. If, for a moment, 


we substitute for _one of these constants, say fo, a variable parameter wu and 
call the new functions simply /,1(0, s, «) and s (0, w) we may write for (51) 


| -s964) | 
Een yr in As So see 
(54). 7 S lim lim > 89 (9, %) bb ode RA a Ou . 
2 (9) s+so(0) wn’ fy ' (0, 8, u) 2wo(6)| O | : 
aul [9, So (0), u| Uh ho 


The partial derivative in the denominator is simply q_'@ (thanks to the choice 

of u=ho, more general we could have chosen w=ho and obtained the same 
tee es l 

1— 6 ae. ji b — 6 AIT ee thi 


in comparing with (13) and (14) proves (16). 


simplicity) and the numerator is 3h, 0 
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Appendix B 


Numerical calculation of the constants h' for M>1. For simplicity we introduce 
constants c’ deriving from c,(y) in the same way as h', from hp (y) (12), thus 


(55) hn=1—en, hi = —e), (121). 


Further we introduce two new constants 


(56) y=VM*—-1, 6 = log Fs 
From (5) we then get 
d' MX 2 
Pa - 2 (M +1) i 
5 é es : 
(57) é Ee [e "Cn (n)| 2M : Pr (14) 


{206 | |/ Mh |! Vena ibe 


Using the Leibnitz formula and putting 7=1 we obtain 


wm S()eer 2 frsnfo (0) a(ean 


Here for abbreviation two functions D and EF have been introduced, having 
the simple form 


(59) D(t)= log (Vi+#2+¢) ; 
14+2¢?? 


- This easily follows frow (5), (6), (56) and (57). 
(58) is a general formula for a successive determination of c!, suitable for nu- 


| merical work. In fact, the integrals in (58) are easily calculated by develop- 
| ing D and E 


Bi ae {27+1 
b>. ; 
; (61) AY) > v Vere 
ae ae 
(62) E()=2t+ > (‘)e y+1) i” 
v=0 
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2 
Thus tD(t) and E(t)—2t¢ are functions of (? and in (58) ? = = has even for 


the smallest M-value in question M=2 a small value <1/3. First we have 
to calculate the terms of the series development of D'(t) B(t) = 2tD'(t) + 
+[E()—2t]D'(t), written as the sum of an even and an odd function. But 
even for the application to M/=2 it is only necessary to compute the first few 
terms of these developments because of the rapid convergence within the inter- 
val of integration in (58). Besides, for increasing values of 7 (m small) the 
integral becomes very small compared to the terms of the left member of (58) 
so that in any case a calculation of the developments up to 7<8 is sufficient. 
In substituting these developments in (58) we get integrals of the type 


1 


+ 
(63) eS i) P,(u) wd (k positive integer). 
=A 


The computations are facilitated by the facts that Zn, vanishes except when 
k—n is a positive even integer or zero, and then we get from the general 
formula of Inx! 


2 y k—n 


(64) Tor = mae hi= 49? and In+i2%=Ink 


(n=0,1,2,...)§ 


The division of the function D'# into the sum of an even and an odd func- 
tion made above is now useful: in fact, when i—n is odd only the part 


2" pi (“ contributes to (58) and when 7—n even only the other part. The 
$a ve 
resulting development in 1/y*® for the integral begins in the former case with 


a term in 1/y*1 if ++12mn and 1/y" if 1+1<n except of course for 7=0 
when the integral vanishes for n= 3, i.e. we have cg=cs;=c,...=0. Fori—n 
even the largest term of the integral is of the order 1/y' if i=n and 1/y" if 
isn. As 6 and 1/y are of the order 1/M it follows from (58) that c’, is of 
the order 1/M' for i=n, and 1/M” for i<n (with the exception above). 
From (21), (22), (24) and (25) it then follows that ani, dni, bi, B; and « all are 
of the order 1/M'. Thus for i<n c’, loses importance, particularly as for in- 
creasing (63) also decreases. This justifies the approximation c,(7)=0 for 
sufficiently large m and it is easy to study in detail its validity in each 
special case. 


Summary 


The stationary space energy distribution of neutrons in an infinite, homo- 
geneous slowing-down medium, containing nuclei of mass M (in units of the 
neutron mass), has been studied. The neutrons are supposed to be emitted 
isotropically with constant energy Ep from a point source. The mean free path A 


' Cf. for instance WHITTAKER and Warson, “Modern Analysis” § 15, 211. 
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is assumed to be constant for energies below Hy. For neutrons of energy EL 
at the distance R from the source (putting A=1) the number of collisions per 
unit volume, per unit energy range is 


—vR 
be’, 


F,(R, E) = — 


where Q is the number of neutrons emitted per unit time and x= log Hy/EH. 
y and to the first approximation also w are functions of R/x only. Analytic 
expressions have been calculated for these functions. When R/x varies from 
0 to co, » increases monotonically from 0 to 1, at first for very small R/« 
proportionally with R/x, for intermediate values of R/x more slowly. In parti- 
cular, for R/x = M/2 (for large M-values; larger than this value for small M), 
roughly corresponding to 20.4, in dividing these values of R/zx, up to co, into 
ez + a where a, } and ¢ are constant 
within each of these intervals (a close to 1). jw is approximately constant for 
very small R/z, and then decreasing with increasing R/zx, for intermediate 
values of R/x, roughly as z/R, and, for very large R/z, as (a/R). The for- 
mula is valid at any distance R for neutrons considerably slowed down (#>1), 
but for increasing FR also for higher neutron energies (Rv>1). The assumption 
of constant A is a good approximation in general only for Hy) = 0.1 Mev (roughly) 
but by chosing a suitable mean value for 2 the formula may be applied even 
for higher #,, at least qualitatively. Applications have been made to deuterium 
and carbon. 


some few intervals we have y=a—2|/ 


The problem of the space energy distribution of neutrons has been suggested 
to me by Professor Ivar WALLER, and I wish to express my gratitude to him 
for continued interest and encouragement and for valuable criticism during 
the work. 
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and Mr. Lars Rupstro6m for assisting me in carrying out this numerical work. 
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